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This paper presents a measurement of the charged current interaction rate of the electron neutrino 
beam component of the beam above 1.5 GeV using the large fiducial mass of the T2K n° detector. 
The predominant portion of the v e flux (~85%) at these energies comes from kaon decays. The 
measured ratio of the observed beam interaction rate to the predicted rate in the detector with 
water targets filled is 0.89 ± 0.08 (stat.) ± 0.11 (sys.), and with the water targets emptied is 
0.90 ± 0.09 (stat.) ± 0.13 (sys.). The ratio obtained for the interactions on water only from an 
event subtraction method is 0.87 ± 0.33 (stat.) ± 0.21 (sys.). This is the first measurement of the 
interaction rate of electron neutrinos on water, which is particularly of interest to experiments with 
water Cherenkov detectors. 


I. INTRODUCTION 

This paper reports a measurement of the ratio of the 
charged current v e event rate relative to the simulation 
with NEUT PQ event generator, version 4.1.4.2, for neu¬ 
trino energies above 1.5 GeV in the T2K beam. The 
interaction rate of electron neutrinos on water has never 
been measured at the neutrinos energies above 1.5 GeV, 
or at any other energies. The mean reconstructed energy 
of the selected neutrinos in the analysis presented in this 
paper is 2.7 GeV. The v e cross section has been measured 
on a liquid freon target for energies between 1.5 GeV and 
8 GeV by Gargamelle [2] and on 12 C for energies around 
32 MeV at LANSCE [3] . Also at lower energies, the 
anti-electron neutrino interactions have been measured 
by experiments near nuclear reactors. A review of neu¬ 
trino cross section measurements can be found in [4|. 

The T2K experiment [5] was built with the primary 
goals of precisely determining the oscillation parameter 
013 via electron neutrino appearance, and of the pa¬ 
rameters 023 and Ato§ 2 via muon neutrino disappear¬ 
ance. The predominantly beam for these measure¬ 
ments is produced at the Japan Proton Accelerator Re¬ 
search Complex (J-PARC) in Tokai. The neutrinos from 
this beam are observed at a near detector, ND280, which 
is located 280 nr downstream from the production tar¬ 
get, where the neutrinos are not expected to have been 
affected by oscillations. The T2K far detector, Super- 
Kamiokande (SK), then measures the muon and electron 
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neutrinos (and anti-neutrinos) after they have undergone 
a near maximal oscillation. 

The oscillation probability for v M i v e depends on the 
mixing parameter, 0 i 3 , and on sub-leading effects that 
depend on the CP-violating phase, <5 cp > and on the mass 
hierarchy [B]. T2K has already observed the appearance 
of 28 v e candidate events at the far detector with a 7.3 a 
significance over a background expectation of 4.92 ± 0.55 
events for 0 13 = 0 jTj. The largest irreducible back¬ 
ground for the appearance measurement comes from the 
predicted 3.2 intrinsic v e beam events. 

In T2K the v e are expected to represent about 1.2% of 
the total neutrino flux [5] . The T2K beam is produced 
by magnetic focusing of pions and kaons produced by 
the interaction of a proton beam with a graphite target. 
The unavoidable i/ e component comes from the decay of 
muons from pion decay, and from kaon decay. In any 
long-baseline neutrino experiment proposed to measure 
CP violation and precisely measure neutrino oscillation 
parameters, the v e component of the beam will be the 
main source of background EHEEj. 

The measurement of the beam v e charged current 
(CCz/ e ) interactions on a plastic scintillator and water 
target using ND280 tracker, was reported in[12]. This 
paper reports a direct measurement of this component 
of the charged current (CC) neutrino interactions in the 
ND280 7r° Detector (P0D) [E], which is located just 
upstream of the tracker. In this selection, the majority 
of the electron neutrinos were produced in kaon decay, 
and have energies above 1.5 GeV. The P0D detector has 
water targets that can be filled or emptied. Data were 
taken both with the targets filled to create a water tar¬ 
get (water configuration), and empty to leave just air in 
place of the water target (air configuration). With data 
in the two configurations a subtraction analysis obtained 
the interaction rate just on water. 

Similar to the subtraction analysis presented here, a 
ratio analysis has been conducted by the Minerva col¬ 
laboration for 2-20 GeV on C, Fe, and Pb compared 
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to CH[T3]. A subtraction analysis of the Minerva data 
is presented the thesis of B.G. Tice [T5]. Apart from 
the Minerva measurements, this appears to be the only 
other use of the subtraction analysis to date in neutrino 
scattering experiments. 

The v e and come from the same pion to muon 
to electron decay chain, and lepton universality allows 
the expected rate of v e to be constrained by measuring 
the much larger flux of v Details concerning the T2K 
beam flux measurement, and further information on re¬ 
cent measurements of v„ interactions in the near detec¬ 
tors, can be found in Ref. [Hi- 

One of the systematic uncertainties in long base¬ 
line neutrino oscillation measurements using water 
Cherenkov detectors comes from model uncertainties in 
the meson exchange current for C versus for O. Having 
measurements of neutrino interaction rates on water is 
therefore important. For a recent review of cross sec¬ 
tion measurements on various nuclear targets refer to the 
PDG[T7]. The only measurements of v tl neutrino interac¬ 
tions on water were reported by the K2K experiment for 
quasi-elastic interactions |18l , and for reactions resulting 
in pions in the final state jTW422j . 

The paper is organized as follows. In Section [TT] the 
P0D detector, used to do the measurement is described. 
The electron selection, and expected backgrounds are 
then described in Section m The particle identification 
(PID) to select electrons from muons in the P0D is a key 
component of this measurement, and will be described 
further in the section on event selection. The water sub¬ 
traction method is then described in |IV| The detector, 
reconstruction, flux and cross section systematic uncer¬ 
tainties in the measurement are reviewed in Section El 
Finally the results of the rate measurement are presented 
in Section |VI| and a summary is in Section Em 

II. ND280 7 r° DETECTOR 

The T2K ND280 7r° Detector (P0D) is a scintillator 
based tracking calorimeter optimized to measure neutral 
current 7r° in the momentum range that contributes to 
backgrounds for v e appearance m ■ Refer to Fig. 4 of 
[T3| for a picture of the P0D detector. The P0D is com¬ 
posed of layers of plastic scintillator alternating with wa¬ 
ter targets and brass sheets or lead sheets. The P0D sits 
in front of a tracking detector made up of two fine grain 
scintillator modules which serve as active targets placed 
between three time projection chambers. Both the P0D 
and tracking detector are surrounded by electromagnetic 
calorimeters and are in a 0.2 T magnetic field. 

The P0D is constructed using 40 scintillator modules, 
each module is constructed with two perpendicular ar¬ 
rays of triangular scintillating bars and is approximately 
38 mm thick. The scintillator modules are arranged in 
three regions. The most upstream and downstream re¬ 
gions of the detector are composed of seven modules in¬ 
terleaved with 4.5 mm thick sheets of stainless steel-clad 


TABLE I. Summary of T2K runs and the number of protons 
on target (POT) used in the analysis. 


T2K run 

P0D Config. 

Beam Power (kW) 

POT (xlO 19 ) 

Run I 

Water 

50 

2.96 

Run II 

Water 

120 

6.96 

Run II 

Air 

120 

3.59 

Run III 

Air 

178 

13.5 

Run IV 

Water 

178 

16.5 

Run IV 

Air 

178 

17.8 

Total 

Water 


26.4 


Air 


34.9 


lead that function as 4.9 radiation length electromagnetic 
calorimeters to improve the containment of photons and 
electrons. The central region serves as a target contain¬ 
ing water. It has 25 water target layers that are 28 mm 
thick sandwiched between 26 scintillator modules and 
1.3 mm brass sheets, positioned in between water targets 
and scintillator layers. The target region has a fiducial 
mass of approximately 1900 kg of water and 3570 kg of 
other materials. 

The energy resolution of the P0D can be estimated 
from Monte Carlo studies by calculating the difference 
between true and reconstructed energy for many events. 
The energy resolution for electrons after the selections 
described in |III| is 16%. 

III. EVENT SELECTION 

A. Overview 

In this analysis, all the data collected between January 
2010 and May 2013 except for very small fraction of Run 
III data, due to the magnetic horn current decrease which 
caused a failure in good spill pre-selection, are used. The 
data are subdivided into different run periods and P0D 
configurations as shown in Table [IJ The simulated data 
used in this analysis corresponds to ten times the Protons 
on Target (POT) of the data, and reproduces the vari¬ 
ous experimental conditions of the different data-taking 
periods. 

Neutrino interactions in ND280 are simulated with the 
NEUT pp event generator, version 5.1.4.2. The genera¬ 
tor covers a range of neutrino energy from several tens 
of MeV to hundreds of TeV and simulates all the nu¬ 
clear targets present in ND280. In the simulated data, 
neutrino interactions are generated outside and within 
the full ND280 volume including all active and inactive 
material, providing information to understand the signal 
and backgrounds from interactions outside the ND280 
fiducial volume. The details of the simulation process 
are described in [28] . 

Simulation of products of the neutrino interactions in 
the P0D is done using a GEANT 4.9.4 simulation [TT| - 
EZ|. The standard GEANT physics list for electromag- 










5 


netic interactions is used in the simulation. 

The analysis uses two reconstructed objects, a track 
and a shower. Within the P0D reconstruction algorithm, 
hits in P0D scintillator layer associated with a recon¬ 
structed track classified as an electromagnetic track (typ¬ 
ically electrons or photons) are forwarded to the shower 
reconstruction stage. Hits associated with a track that 
are classified as a light track (typically muon) or a heavy 
track (typically proton) are not forwarded to the shower 
reconstruction stage and cannot be reconstructed as a 
shower. 

The signal events for the analysis are the charged cur¬ 
rent v e interactions in the P0D. A cut-based event se¬ 
lection using known reconstruction characteristics was 
tuned to maximize the product of efficiency and purity. 
To avoid bias, the selection strategy was developed based 
on Monte Carlo (MC) samples. Event displays of a typ¬ 
ical CCv e candidate and a 7r° background event selected 
in the analysis are shown in Fig. [l] 


B. Selection Cuts 

The event selection strategy focuses on identifying sin¬ 
gle high-energy electron shower events with a vertex in 
the P0D. As a pre-selection, the reconstructed shower in 
the P0D must be in time with the beam bunch time. The 
P0D reconstruction searches for both tracks and show¬ 
ers with two independent algorithms, and the highest 
energy track and the highest energy shower are used in 
the analysis. The reconstruction algorithm builds tracks 
and showers from hits, but as the shower reconstruction 
occurs after the track reconstruction the algorighm needs 
to make sure that the hits shower reconstruction uses are 
the same hits the track reconstruction uses, for each sin¬ 
gle event. Therefore 80% of the hits associated with the 
track and shower are required to be the same. 

In addition, events are selected where the angle of the 
reconstructed shower with respect to the z-axis, which is 
approximately the beam axis, is less than 45°. The scin¬ 
tillator bars of the P0D have a triangular profile with 
angles of approximately 45°. Particles with an angle of 
more than 45° with respect to the beam axis would there¬ 
fore hit more than two adjacent bars in a layer. The 
P0D reconstruction algorithm currently only handles up 
to two adjacent bar hits in a layer, causing reconstruction 
failures for higher angle tracks. 

For this analysis, only events with a reconstructed neu¬ 
trino energy of 1.5 GeV or more are selected. Recon¬ 
structed neutrino energy is calculated from the recon¬ 
structed electron energy and the electron angle using the 
quasi-elastic approximation. In this energy region, the 
majority of the v e flux arises from kaon decays and the 
P0D shows good performance to distinguish electrons 
from other particles. In addition, using a high neutrino 
energy cut improves the purity of the electron sample. 

To reject muons, the median width of the selected 
track is used. In each scintillator layer, the energy- 
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FIG. 1. Side view of a CCu e event (top) and a tt° background 
event (bottom) reconstructed in the P0D. Triangles are hits 
colored by the charge deposited, the green cross symbol shows 
the reconstructed shower vertex, and the green dashed lines 
show the cones of reconstructed showers. 


weighted standard deviation of the position of the hits 
reconstructed in the track is calculated as follows: 

1. If the two hits with the highest deposited energy are 
in adjacent strips, replace them with a single hit. 
The new hit’s position is at the energy-weighted 
average position of the two original hits, and its en¬ 
ergy is the sum of the energies of the original hits. 
Any other hits in the layers are left unchanged. 
This procedure gives layers with minimum ionizing 
tracks very small (almost always zero) width. 

2. The energy-weighted standard deviation of the hit 
positions is calculated for each layer. 
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Track Median Width (mm) 



Track Median Width (mm) 


FIG. 2. Distribution of events which pass all the selection 
criteria with the exception of the track median width cut, for 
water (top) and air configuration (bottom). The vertical line 
shows the cut value used (1 mm). A sudden drop of events 
above 11 mm is an effect of shower median width cut. 


3. Median width is the width of the middle layer after 
ordering by layer width. 

The design of the P0D with layers of high density ma¬ 
terials (brass and lead) causes electrons to shower. The 
reconstructed track of an electron is therefore typically 
wider than the reconstructed track of a muon. This fea¬ 
ture can be used to distinguish muons and electrons with 
the median width of the reconstructed candidate track. 

The track median width for events which pass all the 
selection criteria with the exception of the track median 
width cut, is shown in Fig. [2] and indicates that most of 
the background muon events are rejected by this cut. 

Similarly, to reject background events that contain 
neutral pions, a cut is applied to the median width of 
the selected shower. The shower reconstruction looks for 
hits in a cone from the reconstructed vertex position and 
combines them in one or more showers. It can happen 
that hits from several particles are combined in one recon¬ 
structed shower, especially when they are almost overlap¬ 



Shower Median Width (mm) 



Shower Median Width (mm) 


FIG. 3. Distribution of events which pass all the selection 
criteria with the exception of the shower median width cut, 
for water (top) and air configuration (bottom). The verti¬ 
cal line shows the applied cuts which are optimized for each 
configuration. 


ping. The P0D v e analysis looks for events with a single 
electron. Events with a very wide candidate shower are 
rejected, because such events are more likely background 
events with several particles. The shower median width 
is calculated the same way as the track median width. 
Distributions of events which pass all the selection crite¬ 
ria with the exception of the shower median width cut 
is shown in Fig. [3j It shows many 7r° background events 
are rejected with this cut. 

Finally, a cut is applied to the fraction of the event’s 
charge that is contained in the selected shower. To select 
CCv e events with a high purity, the fraction of the event’s 
charge contained in the candidate shower is exactly 1.0 is 
required, which selects only events with a single shower 
and without muon-like tracks in final state. 
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TABLE II. The selected number of MC signal events, MC 
background events, and the total number of selected MC 
events normalized to data POT for water and air configu¬ 
ration are listed together with the selected data events. In 
addition, the water configuration MC events are split up in 
on-water and not-water events. The errors correspond to the 
statistical uncertainty due to the limited MC statistics. 



MC Signal MC Background 

MC Total 

Data 

Water 

196.1 ±4.8 

56.7 ±2.7 

252.8 ± 5.5 

230 

On-Water 

60.2 ±2.6 

14.5 ± 1.3 

74.7 ±2.9 


Not-Water 

135.9 ±4.0 

42.2 ±2.3 

178.2 ± 4.6 


Air 

173.6 ±4.6 

97.4 ±3.6 

271.0 ±5.8 

257 


TABLE III. The signal efficiencies e and purities p are listed 
for water and air configuration. Events of the POD water con¬ 
figuration are split into events happening on-water and not- 
water. The errors correspond to the statistical uncertainty 
due to the limited MC statistics. 



Efficiency e 

Purity p 

Water 

(10.9 ±0.3)% 

(77.6 ±2.5)% 

On-Water 

(9.8 ± 0.4)% 

(80.6 ± 4.7)% 

Not-Water 

(11.5 ±0.4)% 

(76.3 ±3.0)% 

Air 

(11.0 ±0.3)% 

(64.1 ±2.2)% 


C. Selected Event Samples 

The selected number of events passing all cuts pre¬ 
dicted by the simulation, both when the P0D is config¬ 
ured to contain water and air, together with the number 
of selected data events are presented in Table [TT] The 
water configuration simulation events are separated into 
on-water and not-water events. On-water events are de¬ 
fined as events with true interaction vertex in the water, 
and not-water events have the true interaction vertex on 
scintillator, lead, brass, or other materials besides wa¬ 
ter. All events in the air configuration MC are not-water 
events as the water targets are drained. 


D. Efficiency and purity 



FIG. 4. Selection efficiency of signal events as function of the 
true neutrino energy E tT ue for water and air configuration. 
The error bars correspond to the uncertainties due to limited 
MC statistics. 


IV. WATER SUBTRACTION METHOD 

The measured v e interactions that were collected dur¬ 
ing POD water and air configuration running are com¬ 
pared with the number of v e interactions predicted by 
the POD water and air configuration MC, respectively. 
The measured number of v e interactions are extracted by 
subtracting the predicted MC background B from the se¬ 
lected data events D, resulting in: 


N, 


Data 


= i). 


- ik 


and 


-\jData 

-^CCz^air 


-^air -^air • 


(1) 

( 2 ) 


The background subtracted data are then divided 
by the predicted Monte Carlo signal S to obtain the 
data/MC ratios for the water and air configurations: 


R* 


N uaza 
^'CCt'e ,wate 


, and 


(3) 


R, 



(4) 


The efficiency e and purity p of the simulated electron 
neutrino signal events, for water and air configurations, 
are summarized in Table III In the POD water config¬ 
uration, events are split into events happening on wa¬ 
ter (on-water) and events on scintillator, brass, and lead 
(not-water). 

The selection efficiency of signal events as function of 
the true neutrino energy E tlue for POD water and air 
configurations are shown in Fig. [4] The selection of low 
energy signal events is suppressed by the high neutrino 
energy cut at 1.5 GeV while the selection of high energy 
signal events is suppressed by the shower median width 
cut and the shower charge fraction cut. 


To extract the measured number of on-water charged 
current v e interactions, the measured CC v e interactions 
with POD water and air configurations are compared by 
taking into account the different collected POT and the 
different reconstruction efficiencies for the water and the 
air data sample using: 


j\t Data 

iV COe ,on-water 


( D -water B water ) 

^not-water ’ POT wa ter 

e air • POT air 


(A* 



























In this formula, POT water = 2.64 x 10 20 (POT air = 
3.49 x 10 20 ) is the collected data POT for the P0D water 
(air) configuration. The resulting data/MC ratio for on- 
water CCv e interactions is given by: 


R( 


•on-water — 


AjData 

v CCp e ,on-water 
^on-water 


( 6 ) 


V. SYSTEMATIC UNCERTAINTIES 


The systematic uncertainties in the measurements are 
divided into three categories: detector, reconstruction, 
and neutrino flux/cross section uncertainties. Control 
sample events to study systematic effects in the measure¬ 
ment have been studied, but often the events in these con¬ 
trol samples are not used for the final systematic uncer¬ 
tainty evaluation. The control sample events were found 
to be too similar to the signal events, or did not have the 
same background as the signal events. For this reason a 
simple KS test is used for several of the systematic uncer¬ 
tainty tests, particularly where no deviation is indicated 
in the test. 


A. Detector systematic uncertainties 


The detector’s as-built mass and its mass in the Monte 
Carlo are different. The masses for water and air config¬ 
urations as well as different run periods also vary. These 
differences are incorporated in the analysis procedure 
by re-weighting MC events with mass uncertainties es¬ 
timated to be 0.01 for all configurations. Similarly, the 
fiducial volume and the alignment of the P0D is con¬ 
sidered. Varying the fiducial volume by the MC vertex 
resolution and shifting in P0D alignment provides an es¬ 
timate of the systematic uncertainties in data/MC ratios. 
The uncertainties obtained are smaller than 0.01 for all 
ratios making them negligible in this measurement. 

Possible systematic effects on the reconstructed elec¬ 
tron energy are also studied. The effects are investigated 
by changing the reconstructed energy scale to observe the 
differences in CCv e data/MC ratios. The possible effects 
are as follows: 1. P0D material density and thickness, 2. 
drifts in the P0D response over time, and 3. the simula¬ 
tion (GEANT4) uncertainty in the electron energy depo¬ 
sition. It is assumed the water and air configuration are 
correlated for the P0D material density and thickness 
only. The resulting systematic uncertainties for water 
( d’water )• air ( V a j r ), and on-water (.Ron-water) are 0.05, 
0.05, and 0.10 respectively. 


B. Reconstruction systematic uncertainties 

1. Track PID 


As described earlier at the beginning of Section III the 
classification of the reconstructed tracks is based on the 
P0D PID. Differences in the PID between data and MC 
can therefore cause systematic uncertainties in the CC^ e 
data/MC ratios. 

A PID study with stopping muons in the P0D was 
performed to estimate this uncertainty, and a map of 
mis-PID between a data sample and a simulation of stop¬ 
ping muons was constructed. To estimate the impact of 
the track PID uncertainty on the CCv e data/MC ratios, 
the MC signal and background was weighted according 
to the uncertainty of the map. The systematic parame¬ 
ter values were randomly varied assuming that the water 
and air samples are uncorrelated and also that the sig¬ 
nal and background uncertainties are uncorrelated. The 
uncertainties for water (iZ W ater)j air (I? a „■), and on-water 
(Ron -water) were determined to be 0.05, 0.05, and 0.09 
respectively. 


2. Track and Shower Median width 

To estimate the systematic uncertainty caused by the 
track median width, the plots with all selection criteria 
applied but failing the track median width cut (the N-l 
plots) are integrated, and a Kolmogorov-Smirnov test is 
performed to test if the data and the Monte Carlo event 
distributions are consistent [28] [29]. The Kolmogorov- 
Smirnov test returns a p-value of 91.2% for water and 
92.2% for air configuration indicating that there are no 
significant evidence for a shift between the data and MC 
event distributions. The systematic uncertainty due to 
the track median width cut is therefore negligible for this 
analysis. 

The threshold of the shower median width cut is placed 
in a region with a large number of events. The system¬ 
atic uncertainty on the measured shower median width 
therefore has a larger impact on the CC^ e data/MC 
ratios than the track median width uncertainty does. 
To estimate the systematic uncertainty caused by the 
shower median width, the N-l plots are integrated, 
and a Kolmogorov-Smirnov test is performed. The 
Kolmogorov-Smirnov test returns a p-value of 50.0% for 
water and 65.9% for air configuration. To determine a 
reasonable scaling factor range for Monte Carlo, differ¬ 
ent scaling factors from 0.9 to 1.1 were applied to Monte 
Carlo and the resulting p-values were studied. For a p- 
value of 68%, the peak scaling factor ranged from 0.98 to 
1.02. The systematic effect on the CCu e data/MC ratios 
for -Rwater, Rair, and R on - water coming from the shower 
median width are estimated by varying the scaling fac¬ 
tor that is applied to the MC shower median width. The 
uncertainties obtained for R wa ter, Rain and -R 0 n-water are 
0.04, 0.04, and 0.08 respectively. 
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3. Shower Charge Fraction 

To estimate the possible impact of systematic effects 
of the shower charge fraction on the analysis, additional 
reconstructed objects with low energy are studied. Such 
additional tracks or showers would cause an event to fail 
the shower charge fraction selection criteria. Looking at 
the event distribution of these events, the only hint for 
a systematic difference between data and MC appears in 
the highest bin of the air configuration. Events with a 
shower charge fraction between 0.98 and 1.00 which pass 
all other selection criteria are analyzed to estimate the 
systematic uncertainty. The data/MC difference in this 
region is considered to be the uncertainty on the MC 
events in the signal region, resulting in the systematic 
uncertainties for i? water , -Rain and R 0 n _ W ater of 0.01, 0.04, 
and 0.04 respectively. 


TABLE IV. Summary of systematic uncertainties on the CC^ e 
data/MC ratios for water (R„ a ter), air (f? a n), and on-water 


-on-water 


Systematic Uncertainty 

-Rwater 

-Rair 

Ron-water 

MC Statistics 

0.03 

0.04 

0.12 

P0D Mass 

0.01 

0.01 

0.01 

P0D Fiducial Volume 

< 0.01 

< 0.01 

< 0.01 

P0D Alignment 

< 0.01 

< 0.01 

< 0.01 

Energy Scale 

0.05 

0.05 

0.10 

Hit Matching 

< 0.01 

< 0.01 

< 0.01 

Track PID 

0.05 

0.05 

0.09 

Energy Resolution 

< 0.01 

< 0.01 

0.01 

Angular Resolution 

< 0.01 

< 0.01 

0.01 

Track Median Width 

< 0.01 

< 0.01 

< 0.01 

Shower Median Width 

0.04 

0.04 

0.08 

Shower Charge Fraction 

0.01 

0.04 

0.04 

Flux and Cross Sections 

0.07 

0.09 

0.06 

Total 

0.11 

0.13 

0.21 


C. Flux and cross section systematic uncertainties 


For the inclusion of the flux and cross section system¬ 
atic uncertainties in the analysis, each analyzed MC event 
is re-weighted according to the uncertainties of the flux 
and cross section parameters which are correlated. The 
parameter values and uncertainties are provided by dif¬ 
ferent external measurements such as NA61 and other 
hadronic production experiments, and these parameters 
are then fitted to ND280 data from TPC and FGD, the 
other subdetectors of ND280 than P0D. The systematic 
parameters and their uncertainties obtained from the fit 
to the ND280 data, which includes 25 flux parameters, 
6 FSI parameters, 2 NEUT parameters, and 13 neutrino 
interaction parameters, has been studied in Ref. 1231 . 

To obtain the flux and cross section systematic uncer¬ 
tainties, the systematic parameters are thrown according 
to the covariance matrix and the analysis described in 


Section [TV] is then applied to each throw. The distribu¬ 
tions are fit with single Gaussians and the resulting width 
is considered to be the flux and cross section systematic 
uncertainty for the analysis. The uncertainties obtained 
for water (-R wa ter), air (i? a ir), and on-water (f? Q n -water) 
are 0.07, 0.09, and 0.06 respectively. 


D. Summary of the systematic uncertainties 

All systematic uncertainties on the CC v e data/MC 
ratios for water (i? W ater)> air (i? a - a ), and on-water 
(-Ron-water) that were estimated in the previous sections 
are summarized in Table ED This table also shows the 
total systematic uncertainty. 


VI. RESULTS 

The results obtained for the background subtracted 
data/MC ratio ( R ) for water configuration, air configu¬ 


ration, and on-water are: 


-Rwater = 0.89 ± 0.08 (stat.) ± 0.11 (sys.), (7) 

-Rair = 0.90 ± 0.09 (stat.) ± 0.13 (sys.), and (8) 
Ron-water = 0.87 ± 0.33 (stat.) ± 0.21 (sys.). (9) 

The ratios are consistent with 1, within statistical and 
systematic uncertainties. For the on-water ratio, uncer¬ 
tainties are relatively large due to limited statistics and 
the impact of the subtraction method. 

For the selected events, the distribution of the recon¬ 
structed particle directions is shown in Fig. [5] and the 
distribution of particle energies is shown in Fig. [6] This 
result indicates that the beam i/ e component in high en¬ 
ergy region measured in the data is consistent with ex¬ 
pectations after including constraints from the ND280 
data for all configurations. 


VII. CONCLUSION 

In conclusion, measurements of CCu e interactions us¬ 
ing the ND280 P0D have been made. The P0D includes 
tillable water targets which allows separate measurements 
for the water and air configurations of the ND280 P0D 
as well as the measurement of v e on-water interactions 
above 1.5 GeV in a predominantly beam. About 
~85% of the selected sample comes from the decay of 
kaons. 

The 230 (257) water configuration (air configuration) 
electron neutrino candidate events selected in the data 
are in good agreement with the prediction for the wa¬ 
ter configuration, the air configuration, and for the on- 
water subtraction samples respectively. The measure¬ 
ment is statistically limited, especially for on-water, but 
it will be improved in the future, since collection of ten 
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FIG. 5. Events passing the event selection as a function of 
the particle direction for water (top) and air configuration 
(bottom). The MC events are normalized to data POT, and 
the fit results from ND280 are applied. 


FIG. 6. Events passing the event selection as a function of 
the particle energy for water (top) and air configuration (bot¬ 
tom). The MC events are normalized to data POT, and the 
fit results from ND280 are applied. 


times more data is planned in the coming years. Fur¬ 
thermore, studies and improvements to the reconstruc¬ 
tion algorithms are being investigated to lower the energy 
threshold, which will lead to the measurement of the v e 
cross section on water. 


This is the first v e interaction rate measurement on 
water in the few GeV energy region. Interactions of v e 
on water are of particular interest for long-baseline neu¬ 
trino oscillation experiments, and atmospheric neutrino 
experiments using water Cherenkov detectors with the 
aim to measure CP violation in the lepton sector. 
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